Revealing Sub-Surface Vibrational Modes by Atom-Resolved 
Damping Force Spectroscopy 
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We propose to use the damping signal of an oscillating cantilever in dynamic atomic force mi- 
croscopy as a noninvasive tool to study the vibrational structure of the substrate. We present 
atomically resolved maps of damping in carbon nanotube peapods, capable of identifying the lo- 
cation and packing of enclosed Dy@Cs2 molecules as well as local excitations of vibrational modes 
inside nanotubes of different diameter. We elucidate the physical origin of damping in a microscopic 
model and provide quantitative interpretation of the observations by calculating the vibrational 
spectrum and damping of Dy@Cs2 inside nanotubes with different diameters using ab initio total 
energy and molecular dynamics calculations. 

PACS numbers: 81.05.Tp, 61.48.De, 68.37.Ps 68.35.Ja 



Detecting surface and subsurface vibrations with 
atomic spatial resolution is a daunting endeavor. Re- 
ported observations of molecular vibrations by Inelas- 
tic Scanning Tunnelling Spectroscopy (IESTS) require 
an electrically conducting substrate [l[. Atomic Force 
Microscopy (AFM) experiments involving ultrasonic os- 
cillations of elastically indented samples [H, 0] can be per- 
formed on electrically insulating systems, but yield sub- 
surface images with nanoscale resolution at best. Build- 
ing on the high spatial resolution and sensitivity of dy- 
namic non-contact AFM 0, [Hj], we introduce Damp- 
ing Force Spectroscopy (DFS) as a non-invasive tool to 
study subsurface structure and vibrational modes in com- 
plex molecular systems at the atomic scale. We have 
chosen carbon nanotube peapods @, 0, Q consisting of 
linear chains of fullerenes enclosed in single-wall car- 
bon nanotubes (SWNTs) [9] as a prominent example 
of supramolecular compounds. Of particular interest 
are (M@C n )@SWNT peapods containing M@C n met- 
allofullerenes, hollow cages of n carbon atoms surround- 
ing the metal atom M^ known for their unusual electronic 
transport behavior [10[. 

Here, we demonstrate that monitoring the damping 
of an oscillating AFM tip provides invaluable informa- 
tion not only about topography, but also the subsurface 
vibrational modes that have not been observed before 
with atomic-scale spatial resolution. Our DFS studies of 
(Dy@C 8 2)@SWNT indicate that the observed damping 
of the tip oscillation depends sensitively on its position 
and host tube diameter of (Dy@C 8 2)@SWNT, in agree- 
ment with extensive molecular dynamics (MD) studies 
reported here. Results of our predictive calculations trace 
back the observed damping to the excitation of local vi- 
brational modes by transferring energy from the oscillat- 
ing AFM tip. This truly mechanical oscillator couples 



to the enclosing nanotube first and subsequently to the 
enclosed molecules, revealing their packing structure. 

(Dy@C82)@SWNT peapods were prepared by en- 
capsulating Dy@Cg2 metallofullerenes in open-ended 
SWNTs at a filling rate of ^60%, as confirmed by 
high-resolution transmission electron microscopy and 
DFS 11]. The peapods were deposited at low cover- 



age onto an insulating flat Si02 surface of a Si substrate 
and observed by a home-built dynamic AFM [ll| , shown 
schematically in Fig. QJa). The AFM, equipped with a 
commercial Si cantilever (spring constant of 34.3 N/m, 
eigenfrequency of 159 kHz) and Si tip (nominal tip radius 
of 20 A), was operated at constant oscillation amplitude 
of 21 — 23 A under ultra-high vacuum (p < lxlO -8 Pa) 
at low temperature (T < 13 K). 

High-resolution AFM topography observations of 
SWNTs and peapods, illustrated in Fig. QJa), are shown 
by contour plots in Fig.[2fa:l-4) [ll|. Our topography re- 
sults indicate that the diameter of the empty SWNT in 
Fig.[2](a:l) is d = 16.2±0.5 A. Short-range interatomic in- 
teractions, probed by the dynamic AFM, provide atomic- 
scale contrast in topography images [12[. The observed 
atomic arran gem ent in this tube identifies its chiral in- 
dex as (18, 5) [13]. Comparing these results to those for a 
peapod of the same diameter, presented in Fig.[2ja:2), we 
find the atomic-scale topography features to be amplified 
with respect to the empty SWNT. This contrast enhance- 
ment can likely be attributed to local stiffening of the 
nanotube wall next to the enclosed Dy@C82 molecules. 

The optimum packing geometry of fullerenes in 
peapods is determined primarily by their diameter (dp) 
and that of the enclosing nanotube (d) M- We expect 
an optimum snug fit at d^dp + 6.4 A [lH. [TBI] , with 
the center-to-center inter- fullerene distance dp-p^dp + 
3.2 A, nearly independent of the nanotube diameter. 
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FIG. 1: (Color) (a) Schematic of an oscillating AFM tip scan- 
ning across nanotube peapods with different diameters d at 
different filling levels, (b) Total charge density of the op- 
timized Dy@Cs2 metallofullerene encapsulated in the (18,0) 
nanotube, superposed with the skeletal structure of the pea- 
pod, (c) Schematic of the coupling mechanism between a 
sharp AFM tip and a flexible substrate. The cantilever stiff- 
ness is represented by the spring connecting the tip and its 
suspension. The elastic substrate is represented by the spring 
connecting a representative atom to its surroundings, (d) 
Mechanism of energy dissipation during the dynamic AFM 
operation in non- contact regime. 
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FIG. 2: (Color) AFM topography and damping in SWNTs 
and peapods of different diameter d. (a) Atomically-resolved 
dynamic AFM observations of the topography (a: 1-4) and 
damping (a:5-8), defined as energy loss per cycle AE and 
represented by the color-coded 3D contour plots. Results of 
simultaneous scans for an empty SWNT with d — 16.2±0.5 A 
are shown in (a: 1,5), those for peapods with d = 16.2±0.3 A 
in (a:2,6), with d = 15.0=b0.5 A in (a:3,7), and with d — 
13.0±0.5 A in (a:4,8). (b) Topographic and damping pro- 
files across the same nanotube peapod. The topographic data 
z(x), represented by the black line, are taken along the black- 
white dotted trace in (a:4). The damping data AE(x), repre- 
sented by the red circles, are taken along the red- white dotted 
trace in (a:8). (c) Relationship between nanotube diameter d 
and the maximum damping energy AEmax- The data points 
are representative for tens of tube samples. The dashed line 
is a guide to the eye. 



With the diameter d F ^8.0 A of Dy@C 82 , a A 
wide SWNT should provide an optimum fit, with the en- 
capsulated molecules separated by ^-^^11.2 A. The di- 
ameter d = 13.0±0.5 A of the narrowest peapod observed 
here, discussed in Fig. [2Ka:4), is smaller than optimum. 
In this case, the enclosed molecules distort locally the 
nanotube wall [l6j], as seen in Fig. [2Ka:4). We find the 
observed undulation period of 11.5±0.5 A to agree with 
the above estimated inter-fullerene distance dp- p. 

More intriguing than the topography are our results for 
the spatial variation of the energy loss AE per cycle of the 
cantilever oscillation, defined as damping and monitored 
in DFS [ll|. The damping data, shown in Fig. [2ja:5- 
8), have been recorded simultaneously with the topogra- 
phy data 0,0]. For better quantitative comparison, we 
present topography and damping data for the same sys- 
tem as line plots in Fig. [2Kb). According to Fig. EKa:5), 
almost no damping is observed in an empty SWNT. In 
comparison to these observations, Fig. [2Ka:6) suggests 
that damping increases significantly in a peapod of the 
same diameter. Decreasing the diameter of the enclosing 
nanotube causes an increasingly snug fit of Dy@Cg2 in 
the peapod. Our data in Figs. [2Ka:6-8) indicate that (i) 
AE increases, (ii) the atomic-scale features within AE 
are amplified, and (in) helical undulations in AE grow 
in magnitude with decreasing diameter of the enclosing 



nanotube. In the case of a snug fit, atomic-scale fea- 
tures in AE show a similar or even larger contrast than 
topography features. Even though both are spatially cor- 
related, we observe the strongest damping slightly off the 
center of the largest topographic protrusions on-top the 
encapsulated fullerenes. 

The topography images in Figs. [2Ka) and EJb) show 
some lateral asymmetry, which may be due to Dy@Cs2 
ordering within the peapod or an artifact [ll|. While 
we have taken precautions to avoid artifacts caused by 
tip-induced nanotube displacement (we scanned in the 
axial direction) and feedback lag (we used a small scan 
velocity), we can not exclude an intrinsically asymmetric 
shape or doping profile of the AFM tip as the cause of 
the asymmetry. Irrespective of this finding, we conclude 
that monitoring the damping of the oscillating tip pro- 
vides information superior to topography imaging and 
allows to discriminate between snug and loose packing of 
fullerenes inside the nanotube. 

In Fig. [2Kc) we plot the maximum energy loss per cy- 
cle AE max observed during a DFS scan as a function of 
the host nanotube diameter d. The data points in the 
diagram are presented for a (13,6) and a (11,8) nan- 
otube with d = 13.0±0.5 A, a (19,0) nanotube with 
d = 15.0±0.5 A, and (20,1) and (18,5) nanotubes with 
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d = 16.2±0.5 A. Our results indicate a universal de- 
pendence of AE max on the diameter d only. From the 
narrowest to the widest nanotube, we find the value of 
AE max to decrease from 71 ±5 meV to 30±6 meV and 
7.9±3.0 meV per cycle of the oscillating tip. 

To uncover the origin of damping in DFS, we per- 
formed ab initio density functional calculations of the 
electronic structure, equilibrium geometry and elastic re- 
sponse of these systems, along with MD calculations of 
the dynamical coupling between the AFM tip and the 
peapod. 

Our total energy calculations indicate that unlike the 
most stable C2 isomer of C82, the near-spherical iso- 
mer of Dy@C82 has a different C2 V symmetry. This 
structure is energetically preferred by at least 0.5 eV 
over any other isomer and likely abounds in our sam- 
ples. The total charge distribution of Dy@Cg2 inside the 
(18,0) nanotube, providing optimum enclosure, is shown 
in Fig.IHb). 

The fundamental concept of an AFM imaging an elas- 
tic substrate is depicted in Fig. [He). The position z t of 
the tip apex atom depends on the position z c of the can- 
tilever mount and the deflection of the elastic cantilever, 
represented by a spring. Similarly, an elastic substrate 
is represented by a different spring, allowing the posi- 
tion z a of the closest surface atom to differ from z s of 
the substrate. The force exerted on the tip by the sub- 
strate depends largely on the closest tip-substrate dis- 
tance h = z t — z a . 

To elucidate the origin of damping in the dynamical 
AFM, we first consider the force F a exerted on a sub- 
strate carbon atom by an approaching AFM tip at height 
h, shown by the solid line in Fig. ffl(d). For a given po- 
sition z c of the cantilever mount, the instantaneous po- 
sition of the tip z t and the substrate atom z a depends 
not only on their mutual distance h, which determines 
the nonlinear interaction force shown by the solid line, 
but also on the compensating forces of the strained can- 
tilever and substrate, shown by the dotted lines, as well 
as the history. In the case of a soft cantilever or substrate, 
as the AFM tip approaches from far away, an instabil- 
ity occurs at point A, causing an abrupt decrease in the 
tip-substrate distance to point B. During the retraction 
cycle, a similar instability occurs at C, causing an abrupt 
increase in the tip-substrate distance to D. We note that 
even the closest-approach point B occurs in the non- 
contact regime, since Kb is beyond the equilibrium tip- 
substrate separation He, characterized by F a = 0. The 
hysteresis corresponding to the shaded area in Fig. Hfd), 
delimited by A,B,C and D, represents the energy dissipa- 
tion in the substrate during an ideal approach-retraction 
cycle. 

To quantitatively analyze, how the energy loss of the 
oscillating tip depends on its position and the tube di- 
ameter, we performed a series of MD calculations [11] . 
We induced an initial perturbation by radially displacing 



a surface atom by 0.3 A, corresponding to the distance 
he — h<E in Fig. QJd), and abruptly releasing it. The 
observed energy loss per cycle AE depends on the fre- 
quency hv of vibrational modes excited by the plucking 
process and the likelihood of exciting them. This energy 
is then dissipated efficiently owing to the high thermal 
conductivity of nanotube systems [17]. Since vibrations 
excited in the cantilever do not change with position, 
spatial contrast in AE max is only due to the substrate. 

In empty SWNTs, we found that plucking excites only 
the radial breathing mode (RBM) in the frequency range 
hv = 17.8—23.1 meV in the systems considered here. The 
energy needed to excite this mode is consistent with AE 
values seen in Fig. [2](a:5). We expect stronger damping 
in peapods due to additional modes introduced by the 
fullerenes. Ultra-soft vibrations of Dy within the C82 
cage at hv = 4.1 — 16.5 meV are caused by the soft 
Dy— C82 interaction and the large mass of Dy. Since 
these modes are decoupled from the remaining modes, 
we will ignore the presence of Dy in the following dis- 
cussion. Among the fullerene modes, librations about 
the center are the softest, followed by axial and off- axis 
motion of the center of mass. Harder fullerene modes in- 
clude the RBM at hv = 29 meV, the quadrupolar defor- 
mation mode near hv = 62.0 meV, and higher multipolar 
modes. Especially in the case of a snug fit in the peapod, 
the nanotube and fullerene modes are strongly mixed. 

For each of the nanotubes studied, we changed the ax- 
ial position of the surface atom plucked from on-top a 
fullerene to in-between fullerenes, and also studied a po- 
sition half-way between the two. Our simulations indi- 
cate that damping is strongest for the last position, due 
to the possibility of exciting librations and axial vibra- 
tions of the fullerene underneath. This off-center posi- 
tion, indicated in the inset of Fig. 09(a), agrees with the 
interpretation of our data in Fig. EJa). 

Our MD simulation results for C82@(18,0) are pre- 
sented in Figs. [31(a) andGJb). The radial trajectory of the 
plucked atom, Ar(t) shown in Fig. 09(a), and that of the 
closest fullerene molecule, ArcM(t) shown in Fig. 0(b), 
suggest that several modes are excited. To understand 
which modes dominate damping, we Fourier analyzed all 
trajectories in our study. The Fourier spectra of Ar(t) 
are shown in Fig. 09(c) and those of ArcM(t) in Fig. 09(d) 
for zigzag nanotubes ranging from (17, 0) to (21, 0). Data 
for intermediate nanotube diameters have been obtained 
by interpolation. 

We consider the narrow C82@(17, 0) peapod to provide 
optimum coupling between the snugly fit fullerene and 
the enclosing nanotube. In this case, AE and its spatial 
modulation should be particularly strong. In contrast to 
an empty SWNT, the vibration spectrum of a surface 
atom in Cg2@(17, 0) is very rich, as seen in Fig. [31(c). 
Besides low- frequency modes with hv < 20.7 meV, it 
contains modes with equally spaced frequencies spanning 
the range < hv < 90.9 meV. In the wider (18,0) nan- 
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FIG. 3: (Color) Calculated dynamics and energy dissipation 
in the Cs2@SWNT system following an initial 'radial pluck- 
ing' at the nanotube surface. Strongest damping occurs when 
plucking an atom located axially ~1 A off the center of the 
nearest Cs2- (a) Radial displacement of the plucked atom as 
a function of time in Cs2@(18,0). (b) Radial displacement of 
the C82 center-of-mass as a function of time in Cs2@(18,0). 
(c) Fourier spectrum of the radial displacement of the initially 
plucked atom as a function of the nanotube diameter d. (d) 
Fourier spectrum of the Cs2 center-of-mass as a function of 
the nanotube diameter d. The insets in (a) and (b) illustrate 
the initial distortion and the quantity monitored. 



otube, the vibration spectrum is similar, but softer. In 
the (19,0) and wider nanotubes, the fullerene-nanotube 
coupling is reduced drastically, as fullerenes may rear- 
range at little energy cost. As seen in Fig.[3^c), decreased 
fullerene-nanotube coupling causes a strong intensity 
drop of the high-frequency modes at hv > 20 meV. Es- 
pecially in the narrow nanotubes, our results for the off- 
center plucking position in Fig. [3tc) differ from those of 
other simulations, where the plucked atom is on top of a 
fullerene. Since plucking on top a fullerene excites neither 
librations nor axial fullerene motion, energy dissipation 
in this plucking position is lower. Comparing the results 
for the two plucking positions explains the high spatial 
modulation of AE seen in Fig. [2](a:8). 

The vibration spectrum of the center-of-mass motion 
of a C$2 fullerene in the (17,0) nanotube, shown in 
Fig. [3^d), is rather featureless and much softer than that 
in Fig. [3](c). The hv~8 meV off-axis vibration mode, 
which dominates the (17, 0) spectrum, red shifts in wider 
nanotubes. Owing to the zigzag or helical arrangement 
of fullerenes in wider peapods, the off-axis vibrations of 
the fullerenes may couple to soft axial vibrations and li- 
brations, which appear in the Fourier spectrum of Atcm 
in Fig. [31(d). These important soft modes, which do not 
change Atcm-, modify the spectrum indirectly by cou- 
pling anharmonically to off-axis modes. Especially in the 



wider nanotubes, we find the center-of-mass vibrations of 
C82 to be rather insensitive to the axial position of the 
plucked atom, thus contributing very little to the spatial 
modulation of AE. 

In conclusion, Damping Force Spectroscopy can reveal 
geometrical packing and vibrational spectra of subsurface 
structures at the atomic scale. We expect that this tech- 
nique will gain significant appeal for atomic-level inves- 
tigations of three-dimensional supramolecular systems. 
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